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A novel information-geometric approach to chaotic dynamics on curved statistical manifolds based 
on Entropic Dynamics (ED) is suggested. Furthermore, an information-geometric analogue of the 
Zurek-Paz quantum chaos criterion is proposed. It is shown that the hyperbolicity of a non- 
maximally symmetric 6A^-dimensional statistical manifold Ms underlying an ED Gaussian model 
describing an arbitrary system of 3A'' non-interacting degrees of freedom leads to linear information- 

- - - geometric entropy growth and to exponential divergence of the Jacobi vector field intensity, quantum 

00 ' and classical features of chaos respectively. 
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I. INTRODUCTION 



Entropic Dynamics (ED) [1] , namely the combination of principles of inductive inference (Maximum relative Entropy 
^' Methods, [2]) and Information Geometry (IG) [3], is a theoretical framework constructed on statistical manifolds and 
I < it is developed to investigate the possibility that laws of physics, either classical or quantum, might reflect laws 
1^ ' of inference rather than laws of nature. This paper is a follow up of a series of the author's works [4, 5]. In 
this article, we use the ED theoretical framework to explore the possibility of constructing a unifying (classical 
' and quantum) criterion of chaos. We assume the system under investigation has 3N degrees of freedom, each one 
I— —I: described by two pieces of relevant information, its expectation value and its variance (Gaussian statistical variables). 

This leads to consider an ED model on a non-maximally symmetric 6N- dimensional statistical manifold Ms- The 
" ' manifold Ms has constant negative Ricci curvature proportional to the number of degrees of freedom of the system, 
— — 3A^. An information-geometric analog of the Zurek-Paz quantum chaos criterion is suggested. It is shown 
that the system explores statistical volume elements on Ms at an exponential rate. We define a dynamical information- 
\^ , geometric entropy of the system and we show it increases linearly in time (statistical evolution parameter) and 

Ti^J" ' is proportional to the number of degrees of freedom of the system. The geodesies on Ms are hyperbolic trajectories. 
Using the Jacobi-Levi-Civita (JLC) equation for geodesic spread, it is shown that the Jacobi vector field intensity Jm^ 
I ■ diverges exponentially and is proportional to the number of degrees of freedom of the system. Thus, Rmsj '^Ms ^'^d 
00 '. J Ms ^I'S proportional to the number of Gaussian-distributed microstates of the system. This proportionality leads to 
. conclude there exists a substantial link among these information-geometric indicators of chaoticity. 



II. THE ED GAUSSIAN MODEL 

Given two probability distributions, how can one define a notion of "distance" between them? The answer to 
this question is provided by IG. As it is shown in [6] and [7], the notion of distance between dissimilar probability 
distributions is quantified by the Fisher-Rao information metric tensor. We consider an ED model whose microstates 

span a 3iV-dimensional space labelled by the variables |-^| — x'^'^\...., x^^^} with af'^"-' = {^^^\ x"^\ x"^^^, 

a = 1,...., N and xi"'' e M with a = 1, 2, 3 . We assume the only testable information pertaining to the quantities 

xi"-* consists of the expectation values (x'a'^\ and the variance Axi"-* = \l ii^^a'^ ^ V The set of these 




expected values define the 6A^-dimensional space of macrostates of the system. A measure of distinguishability among 
the states of the ED model is achieved by assigning a probability distribution P [X ) to each macrostate 9 where 
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|o| = |'"'^-*^a"'',*-^-'^i"'*| with a = 1, 2,...., A'' and a = 1,2,3. The process of assigning a probabihty distribution to 

each state provides Ms with a metric structure. Specifically, the Fisher-Rao information metric defined in (7) is a 
measure of distinguishability among macrostates. It assigns an IG to the space of states. 



III. THE STATISTICAL MANIFOLD Ms 



Consider an arbitrary physical system evolving over a 3A''-dimensional space. The variables |-^| = 

{x^^\ x^^\...., label the 3A''-dimensional space of microstates of the system. Each macrostate may be thought 

as a point of a 6A7^-dimensional statistical manifold with coordinates given by the numerical values of the expectations 

2\ 

The available information can be written in the form 



(1)^) = (xi")) and (2)^' = A.xi") ^ 
of the following 6A^ information constraint equations, 

+ 00 



dx, 



(1) 



Ax'^ 



(a) 



+ 00 



where ^^H 



and (2)6»^") = Axi"^ with a = 1, 2,...., iV and a = 1,2,3. The probability distributions Pi 



(a) 



are constrained by the conditions of normalization, 



+ 00 



1. 



(2) 



Information theory identifies the Gaussian distribution as the maximum entropy distribution if only the expectation 
value and the variance arc known. The distribution that best reflects the information contained in the prior distribution 

m updated by the information (^^a""*^ , Axi"-*^ is obtained by maximizing the relative entropy 



(e) =-y d^^xp (x|e)iog 



p 




e) 




m (x 





(3) 



where m{X) is the prior probability distribution. As a working hypothesis, the prior m yX j is set to be uniform since 

we assume the lack of prior available information about the system (postulate of equal a priori probabilities). Upon 
maximizing (3), given the constraints (1) and (2), we obtain 



AT 3 

^(^|Q) = nn^i"^ (4-^ 



(4) 



a=l a=l 



where 



a f^a J "a J 
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(5) 



and ^^^9^^^ = Ma"\ ('^)0^^'^ = cri"-'. The probability distribution (4) encodes the available information concerning the 
system. Note that we have assumed uncoupled constraints among microvariablcs xi"'. In other words, we assumed 
that information about correlations between the microvariables need not to be tracked. This assumption leads to the 



3 



simplified product rule (4). However, coupled constraints would lead to a generalized product rule in (4) and to a metric 
tensor (7) with non-trivial off-diagonal elements (covariance terms). Such generalizations would require more delicate 
analysis. Deviations from Gaussian- type information constraints and the presence of a nonvanishing correlation coeffi- 
cient (r = <('^i~/'^)0^2-^2)) ^ J, jg ^jjg correlation coefficient of the two dependent random variables xiand among 

the random microvariables a;!"' are some of the new topics appearing in a forthcoming paper [8]. 



A. METRIC STRUCTURE OF Ms 



The dimensionless line element ds between P (x|e) andp(x 



Q + d&] is given by, 



where 



dXP X 



e 



d\ogP(X 



e) 5iogp(x e) 



(6) 



(7) 



is the Fisher-Rao metric. Substituting (4) into (7), the metric g^i, on Ms becomes a 6A'' x 6N matrix M made up of 
3N blocks M2x2 with dimension 2x2 given by, 



2x2 







2x(ai"))"'y 
with a = 1, 2,...., A'' and a = 1,2,3. From (7), the "length" element (6) reads 



(8) 



N 3 



a=l a=l 



1 .d^(«)2 + ^_^rf^M2 



2 "'"a 



(9) 



We bring attention to the fact that the metric structure of A4s is an emergent (not fundamental) structure. It arises 
only after assigning a probability distribution P 9^ to each state G. 



B. CURVATURE OF M 



The Ricci scalar curvature R is given by, 

R = g'"'R^., (10) 

where g^'^gvp = f^p so that g^'^ — {g^v) ^ ■ The Ricci tensor i?^^ is given by, 

R^lv = dgV — dijT + r^,^r^^ — rj^^r^^. (ii) 

The Christoffel symbols T^^^, appearing in the Ricci tensor are defined in the standard way, 

= {d^9eu + d.g^, - d,g^,) . (12) 
Using (9) and the definitions given above, we can show that the Ricci scalar curvature becomes 

Rms = -3iV < 0. (13) 

From (13) we conclude that Ms is a GAT-dimensional statistical manifold of constant negative Ricci scalar curvature. A 
detailed analysis on the calculation of Christoffel connection coefficients using the ED formalism for a four-dimensional 
manifold of Gaussians can be found in [5]. Furthermore, it can be shown that A^g is not a pseudosphere (maximally 
symmetric manifold) since its sectional curvature is not constant. Considerations about the negativity of the Ricci 
curvature as a strong criterion of dynamical instability and the necessity of compactness of Ms in "true" chaotic 
dynamical systems will appear in [8]. 
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IV. CANONICAL FORMALISM FOR THE ED-GAUSSIAN MODEL 

At this point, we study the trajectories of the system on A4s- We emphasize ED can be derived from a standard prin- 
ciple of least action (Maupcrtuis- Euler-Lagrange-Jacobi-type) [1, 9]. The geodesic equations for the macrovariables 
of the Gaussian ED model are given by, 

with 1-1=1, 2,...,6N. Observe that the geodesic equations are nonlinear, second order coupled ordinary differential 
equations. 

A. GEODESICS ON Ms 

We seek the explicit form of (14) for the pairs of statistical coordinates f^"^). Substituting the explicit 

expression of the Christoffcl connection coefficients into (14), the geodesic equations for the macrovariables /ii"'* and 
cri"^ associated to the microstate xi"' become, 

dr^ ^(") dr dr dr^ [ dr ) + 2ai") [ dr ) ^'"^ 

with a = l, 2,...., A'' and a = 1, 2, 3. This is a set of coupled ordinary differential equations, whose solutions are 

cosh(2/3i"V) -si„h(2/3i"'r) + ^p^ 

_ p(a) cosh(;3(°)r)-sinh(;3(°)r) 
CTa [T) — Ha / ■ 

cosh (2/3i°' t) -sinh (2/3i°' r ) + \ ^ 1 1 

\ y \ y s(3^°') 



(16) 



The quantities si"^ , cl"^ , /S*^"-* are real integration constants and they can be evaluated once the boundary conditions 
are specified. We are interested in investigating the stability of the trajectories of the ED model considered on M.s- 
It is known [9] that the Riemannian curvature of a manifold is closely connected with the behavior of the geodesies 
on it. If the Riemannian curvature of a manifold is negative, geodesies (initially parallel) rapidly diverge from one 

another. For the sake of simplicity, we assume very special initial conditions: si"^ = A, = A e M"*", Ca"^ = 0, 
Va = 1, 2,...., A'^ and a = 1,2,3 . However, the conclusion we reach can be generalized to more arbitrary initial 
conditions. It is worthwhile noticing that, in our case. Ads is a geodesieally complete manifold since every maximal 
geodesic is well-defined for all temporal parameters r. Therefore, M.s represents a natural setting for global questions 
in this Riemannian geometric framework applied to probability theory and the search for a weak criterion of chaos 
can be carried out [8]. 

V. LINEARITY OF THE INFORMATION-GEOMETRIC DYNAMICAL ENTROPY 

Recall that A^^ is the space of probability distributions P {^X labeled by 6 A'' statistical parameters 6. These 
parameters are the coordinates for the point P, and in these coordinates a volume element dVMs reads. 



dVMs = Vgd'^'e = n n -T^^^dfii^^d^^^^- (17) 



1=1 a=l lak ' j 

The volume of an extended region AVm^ (t; A) of Ms is defined by, 

AVMAr;X)''='l[E I I 7^rfMi"Wi") (18) 
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where /ii"'' (r) and (ri"-* (t) arc given in (16). The quantity that encodes relevant information about the stability of 
neighboring volume elements is the the average volume (r; A), 

T 

Vm. (t; A) ^ {AVm. (t; A)), 1 J AVm, (t'; A) dr' ™ e^^^^ (19) 



This asymptotic regime of diffusive evolution in (19) describes the exponential increase of average volume elements 
on Ms- The exponential instability characteristic of chaos forces the system to rapidly explore large areas (volumes) 
of the statistical manifolds. It is interesting to note that this asymptotic behavior appears also in the conventional 
description of quantum chaos where the entropy increases linearly at a rate determined by the Lyapunov exponents. 
The linear entropy increase as a quantum chaos criterion was introduced by Zurek and Paz [10]. In our information- 
geometric approach a relevant variable that can be useful to study the degree of instability characterizing the ED 
model is the information-geometric entropy quantity defined as, 

Sm/= iim logVMAr; X) . (20) 

r— *oo 

Substituting (18) in (19), equation (20) becomes, 

Sm. ^ 3A^Ar. (21) 



The entropy-like quantity (21) is the asymptotic limit of the natural logarithm of the statistical weight {AVms}t 

defined on Ad^ and it grows linearly in time, a quantum feature of chaos. Indeed, equation (21) may be considered 
the information-geometric analog of the Zurek-Paz chaos criterion. Zurek and Paz considered a chaotic system, a 
single unstable harmonic oscillator characterized by a potential V (x) = (A is the Lyapunov exponent), coupled 

to an external environment. In the reversible classical limit, the von Neumann entropy of such a system increases 
linearly at a rate determined by the Lyapunov exponent, 

Canw^ (Zurek-Paz) ^ Ar. (22) 
In general, the von Neumann entropy 'S'q^'^ntum is given by 

'S'quanrum = (p log^ P) = " X] ^^^.2 (23) 

j 

where the normalized {tr (p) = 1) density operator p is defined as 

P = E^i|V',-)(V',|,?K)=A,|V,). (24) 

3 

Notice that the consideration of 3A^ uncoupled identical unstable harmonic oscillators characterized by potentials 
Vi (x) = {Xi = Xj; i, j = 1, 2,..., 3A^) would simply lead to 

Cntu:^^ (Zurek-Paz) 37VAr. (25) 

The resemblance of equations (21) and (25), either in the form or the content is astonishing. A detailed discussion 
about this result and an additional discussion about a possible connection of (21) to the Kolmogorov-Sinai entropy, 
one of the most powerful indicators of chaos in classical dynamical systems, will appear in [8]. 

VI. EXPONENTIAL DIVERGENCE OF THE JACOBI FIELD INTENSITY ON Ms 



{^'m. ir; mZir'" where 



Finally, we consider the behavior of the one-parameter family of neighboring geodesies Tgms (•^) 

^s7■h^lr^l 

AeM+ 

t^a A) - 2Acosh(2Ar)-sinh(2Ar)+ ' 



(26) 



(s cosh (At) - sinh (At) 

85? 



cosh (2At) - sinh (2At) 
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with a = 1, 2,...., N and a = 1,2,3. The relative geodesic spread on a (non-maximally symmetric) curved manifold 
as Ms is characterized by the Jacobi-Levi-Civita equation, the natural tool to tackle dynamical chaos [11, 12], 



where the Jacobi vector field is defined as, 



dT 



(r; A) 



5\. 



(27) 



(28) 



Equation (27) forms a system of 6A^ coupled ordinary differential equations linear in the components of the deviation 
vector field (28) but nonlinear in derivatives of the metric (7). When the geodesies are neighboring but their relative 
velocity is arbitrary, the corresponding geodesic deviation equation is the so-called generalized Jacobi equation [13]. 
Substituting (26) in (27) and neglecting the exponentially decaying terms in and its derivatives, integration of 
(27) leads to the following asymptotic expression of the Jacobi vector field intensity, 



Jms = \\J\\ = {g,n.J''r)' 



3Ne/ 



(29) 



Further details on the derivation of this result for a four-dimensional statistical manifold are in [5]. We conclude that 
the geodesic spread on Ads is described by means of an exponentially divergent Jacobi vector field intensity j a 
classical feature of chaos. In our approach the quantity Aj, 



Aj '= lim — In 

T— »OOT 



■J Ms 


(r) 


J Ms 


(0) 



would play the role of the conventional Lyapunov exponents. In conclusion, we have shown that, 

Rj^^ = -3N, Sms ~ ^NXt, Jms ~ ^A/'e 



At 



(30) 



(31) 



The Ricci scalar curvature Rms > the information-geometric entropy Sms ^^'^ the Jacobi vector field intensity Jms 

are proportional to the number of Gaussian-distributed microstatcs of the system. This proportionality leads to the 
conclusion that there exists a substantial link among these information-geometric measures of chaoticity, namely 



Rms ~ ~ Jms 



(32) 



Equation (32), together with the information-geometric analog of the Zurck-Paz quantum chaos criterion, equation 
(21), represent the fundamental results of this work. We believe our theoretical modelling scheme may be used to 
describe actual systems where transitions from quantum to classical chaos scenario occur, but this will be argued 
elsewhere [8]. 



VII. FINAL REMARKS 



In conclusion, a Gaussian ED statistical model has been constructed on a 6A''-dimensional statistical manifold Ms- 
The macro-coordinates on the manifold are represented by the expectation values of microvariables associated with 
Gaussian distributions. The geometric structure of Ms was studied. The manifold Ms'^sa curved manifold of constant 
negative Ricci curvature — 3iV . The geodesies of the ED model arc hyperbolic curves on Ms- A study of the stabiUty of 
geodesies on was presented. The notion of statistical volume elements was introduced to investigate the asymptotic 
behavior of a one-parameter family of neighboring volumes J^Vms (•^) — i^Ms •^)}agr+- An information-geometric 
analog of the Zurek-Paz chaos criterion was suggested. It was shown that the behavior of geodesies is characterized 
by exponential instability that leads to chaotic scenarios on the curved statistical manifold. These conclusions are 
supported by a study based on the geodesic deviation equations and on the asymptotic behavior of the Jacobi vector 
field intensity Jm^ on A^^. A Lyapunov exponent analog similar to that appearing in the Riemannian geometric 
approach to chaos [14] was suggested as an indicator of chaoticity. We think this is a relevant result since a rigorous 
relation among curvature, Lyapunov exponents and Kolmogorov-Sinay entropy is still under investigation [15] and 
since there does not exist a well defined unifying characterization of chaos in classical and quantum physics due to 
fundamental differences between the two theories [16]. 
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